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Nanoscale forces play an important role in different scanning probe microscopies, most notably 
atomic force microscopy (AFM). In contrast, in scanning near-field optical microscopy (SNOM) 
a light-induced coupled local optical polarization between tip and sample is typically detected by 
scattering to the far field. Measurements of the optical gradient force associated with that optical 
near-field excitation would offer a novel optical scanning probe modality. Here we provide a gener¬ 
alized theory of optical gradient force nano-imaging and -spectroscopy. We quantify magnitude and 
distance dependence of the optical gradient force and its spectral response. We show that the optical 
gradient force is dispersive for single particle electronic and vibrational resonances, distinct from 
recent claims of its experimental observation. In contrast, the force can be absorptive for collective 
resonances. We provide a guidance for its measurements and distinction from competing processes 
such as thermal expansion. 


Control and enhancement of light-matter interaction 
via radiation pressure and optical gradient forces form 
the basis of cavity optomechanics and optical tweezers 
mm- Fundamentally related, the tip-sample interaction 
in electrostatic, magnetic, or atomic force microscopy 
is based on classical electromagnetic or quantum elec¬ 
trodynamic effects (Van der Waals/Casimir interaction) 
mm- Under light illumination, the induced spatially lo¬ 
calized coupled optical polarization between the tip and 
the sample forms the basis of scanning near-held optical 
microscopy (SNOM) OE], with the near-held signal typ¬ 
ically detected by far-held scattering. However, the cou¬ 
pled optical polarization is expected to also give rise to 
an optical gradient force between the tip and the sample 
as illustrated in Fig. ^ [7HS]. Its mechanical detection 
by means of atomic force microscopy (AFM) techniques 
has been proposed as an promising alternative way for 
scanning probe optical nano-imaging and -spectroscopy 

[SHS]. 

Here, we provide a generalized theoretical description 
of the optical gradient force between a scanning tip and 
sample as an imaging contrast mechanism. The spec¬ 
troscopic response of the optical gradient force is found 
to be dispersive for molecular electronic or vibrational 
and other single particle excitations, in contrast to recent 
experimental claims of gradient force nano-spectroscopy 
[13 H]. We hnd that only collective polariton resonances 
can give rise to absorptive spectral force prohle. The ef¬ 
fect is distinct from the accompanying thermal expansion 
due to optical absorption (Fig. [^), which results in ab¬ 
sorptive resonance spectra in all cases [T5HT7] . While the 
optical gradient force is near the thermal cantilever noise 
limit in a room temperature AFM, the effect should be 
detectable under improved force detection or cryogenic 
conditions as established in cavity opto-mechanics [I]- 

Recent experimental efforts on measuring the optical 
near-held gradient force have not been conclusive. For 
example, a force image of a metal bowtie antenna on res¬ 


onance at A = 1550 nm has been interpreted to result 
from an optical gradient force [18j . In contrast, simi¬ 
lar force images on gold split ring resonators were at¬ 
tributed to thermal expansion |19j . For molecular elec¬ 
tronic resonance spectroscopy an absorptive response has 
been observed and attributed to the optical gradient 
force m I14j . in marked contrast to the predictions in 
this work. In contrast, a nearly instantaneous (few ps) 
response from ultrafast experiments |20j and polariza¬ 
tion dependence in imaging plasmonic particles |21] show 
characteristics of the expected optical gradient force be¬ 
havior. However, both experiments lack spectroscopic 
information which would be desirable for a unique as¬ 
signment. The situation is equally confusing regarding 
theory. Most studies so far have treated both tip and 
sample as point dipoles [13 [H, even at close proximity 
where the dipole approximation fails, or only calculated 
the distance dependence of the force without considering 
its spectral response [T3], leaving key open questions on 
experimental feasibility and distinguishing spectral char- 



FIG. 1. (a) Optical gradient force induced by excitation 

of coupled optical polarization between tip {fit) and sample 
(fis), with resulting force F oc fitfis- (b) The optical excita¬ 
tion is accompanied by absorption, with the resulting thermal 
expansion Ah gives rise to a competing force reaction of the 
AFM tip. 
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acteristics. 

In this work, we numerically calculate the optical gra¬ 
dient force between the tip and sample from a sur¬ 
face integral of the Maxwell stress tensor T^- using fi¬ 
nite element electromagnetic simulation (COMSOL Mul¬ 
tiphysics). The force on the tip is given by 

Fi=Y^J TijUjds ( 1 ) 

j ^ 

with tip surface S, surface normal n, and 

Tij = eo ^EiEj — + ~ ■ 

Since the radius of tip apex r ~ 10 nm is much smaller 
than the laser wavelength A, the near-held tip and sam¬ 
ple interaction under plane wave illumination can be an¬ 
alyzed in the quasistatic approximation by assuming the 
probe to be a polarizable sphere with radius r. The 
resulting held distribution on the sample can be ehec- 
tively reduced to an image sphere of radius r, as shown 
in Fig. While the exact geometry of tip and sample af¬ 
fects both the details of magnitude and spectral response 
of the optical gradient force, especially for strong polari- 
tonic resonances, the limiting case of two hnite spheres 
provides enough general insight into the spectral varia¬ 
tion of the force spectrum and its distance dependence. 

Considering incident polarization along the tip axis, 
and rotational symmetry along the z-axis allows sim¬ 
ulation in reduced dimensions to decrease the compu¬ 
tational complexity. For the held simulation, a dense 
mesh with element size of 0.2 nm or less is applied near 
the spheres, and the distant surrounding is meshed with 
maximum element size of 20 nm. A uniform electric held 
E — 10® V/m along the z direction is applied. This 
electric held corresponds to an average laser intensity of 
Iq = 1.33 mW//im^ as used under typical experiment 
conditions [II1E2]. 

We consider two types of resonant processes of the 
sample: i) single particle molecular electronic and vi¬ 
brational resonances, and ii) collective plasmonic reso¬ 
nances. The results are readily generalizable to samples 
with other types of resonances including phonons, exci- 
tons, and related polaritons. For the case of an electronic 
resonance, a sample consisting of Rhodamine 6G (R6G) 
dye molecules is modeled as a single harmonic oscilla¬ 
tor for its dielectric function e(w) = ei(a;) -I- ie 2 {uj), with 
transition energy of 2.42 eV, and line width of 0.41 eV, 
as a best ht to experimental results. A tungsten tip with 
r = 10 nm is used to ensure a hat spectral response of the 
tip in the relevant visible spectral region. For the molec¬ 
ular vibrational resonance, a gold tip sphere is used to 
interact with a sample sphere (both with r = 10 nm) of 
poly(methyle methacrylate) (PMMA) with characteris¬ 
tic carbonyl resonance (G=0) at 1729 cm“^. Finally, to 
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FIG. 2. (a) Electronic resonant force spectrum between a 

tungsten tip and R6G sample (both of radius r = 10 nm) sub¬ 
ject to an optical field E = 10® V/m, separated by z = 2 nm. 
(b) Dielectric function of R6G e{u}) = ei(a;) -|- ie 2 {u}) mod¬ 
eled as a single harmonic oscillator, (c) Corresponding force 
spectra between a gold sphere and PMMA sample of radii 
r = 10 nm. Dispersive force spectrum at separation distance 
z = 20 nm (multiplied by 20, black) compared to force spec¬ 
trum at z = 2 nm (red), (d) Dielectric function of PMMA 
e{uj) = ei{io) + around the carbonyl resonance at 1729 

cm“^. (e) Simulated PMMA absorption at distances z = 20 
nm (black) and z = 2 nm (blue), (f) The distance depen¬ 
dence of the optical force (red) on resonance shows a complex 
distance scaling, neither following simple dipole-dipole power 
law nor exponential scaling. The slight increase of optical ab¬ 
sorption (blue) is due to local field enhancement in the gap 
region. 


study the case of plasmonic resonance, two silver spheres 
of r = 10 nm are modeled as the tip and sample system. 

Fig. shows the resulting force spectrum between the 
tungsten tip and R6G sample at a distance of z = 2 nm. 
The force spectrum is dispersive, and follows the trend 
of the real part of the dielectric function ei(w) shown 
in Fig. Ih- The relatively small force variation between 
23 fN to 25 fN across the resonance is due to the large 
broadband offset as a result of ei ^ £ 2 , which is a charac¬ 
teristic property of the dielectric function for molecular 
electronic resonances. 

Gorrespondingly, Fig. show the force spectrum of 
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FIG. 3. (a) Optical gradient force spectrum of two silver 

spheres with radius r = 10 nm separated by 2 = 20 nm sub¬ 
ject to external optical field E = 10® V/m. (b) Absorption 
spectrum of the two spheres with symmetric resonance at 3.5 
eV. Inset: real and imaginary part of the normalized polariz¬ 
ability Q = (cAg — l)/(£Ag + 2) for a silver sphere. Force and 
absorption spectra follow the imaginary part of a{uj). 


the Au sphere interacting with the carbonyl vibrational 
resonance of PMMA at separation of z = 20 nm (black). 
The force increases by 20 times in magnitude to a peak 
value of F « 1 fN when decreasing the distance to z = 2 
nm (red). Due to the smaller vibrational dipole moment 
compared to the electronic counter part, the force is much 
weaker than the case of R6G but exhibits a larger rela¬ 
tive spectral variation across the resonance from 0.2 fN 
to 1 fN (for z = 2 nm). Irrespective of distance, and sim¬ 
ilar to the electronic resonance case, the force spectra are 
dispersive and follow the real part of the dielectric func¬ 
tion of PMMA ei(a;) (Fig. [§i) . For comparison. Fig. ^ 
shows the absorption spectra at z = 20 nm and z = 2 nm 
in black and blue, respectively. In contrast to the force 
spectra, the absorption spectra show the well known sym¬ 
metric resonance behavior corresponding to the imagi¬ 
nary part of the dielectric function e 2 (w). Except for 
a change in magnitude, peak position and line shape of 
both optical gradient force and absorption spectra are 
invariant with respect to tip-sample distance. The dis¬ 
tance dependence of the optical force on resonance (red) 
is shown in Fig. [^. Due to increasing multipole con¬ 
tributions with decreasing distance, the distance scaling 
follows neither simple exponential nor power law, and is 
far more shallow than what one would expect based on a 
simple dipole-dipole interaction (ex z^). The magnitude 
of absorption (blue) only increases slightly due to local 
field enhancement with decreasing distance. 

Fig. shows the force spectrum of two silver spheres 
(r = 10 nm) separated by z = 20 nm for an applied 
optical field of F = 10® V/m. The plasmonic reso¬ 
nant force spectrum with a slight asymmetric peaks at 
3.5 eV with F^ax — 30 fN. For comparison. Fig. 
shows the absorption spectrum for the two silver spheres. 
The normalized complex polarizability of a silver sphere 
d = (cAg — l)/(eAg + 2) is shown in the inset. As can be 
seen, both optical gradient force and absorption spectra 


follow the imaginary part of polarizability Im{a){u}). 

Fig.|4}i shows the spectral evolution of the force as a 
function of distance. When the separation of the two 
spheres decreases from z = 20 nm to z = 1 nm, the opti¬ 
cal gradient force increases by a factor of ^ 10® reaching 
F « 10 pN. Correspondingly, the resonance frequency of 
the force spectrum red-shifts due to the plasmonic cou¬ 
pling between the two spheres |23]. In contrast, the mag¬ 
nitude of optical absorption does not change significantly 
with separations as shown in Fig. The distance de¬ 
pendence of the optical gradient force on- (red) and off- 
resonance (blue) is plotted in Fig.|^. Both increase with 
decreasing distance approximately following a power law 
cx z", with n « 2, i.e., again less steep than the point 
dipole z^ dependence. 

The different spectral behaviors of the optical gradient 
force can be rationalized based on the dipole-dipole in¬ 
teraction assuming two independently polarized spheres 
when z > 2r, with resulting force proportional to the 
polarizability, i.e., F oc Re{ata*). For a non-resonant 
tip, the force simplifies to F cx Re{as), describing the 
dispersive force spectra for a tungsten tip and a R6G 
sample in the visible, and a gold tip and a PMMA sam¬ 
ple in the infrared. It is interesting to note that the 
optical gradient force can in principle also become re¬ 
pulsive for Re{as) crossing zero. When tip and sample 
consist of the same material (e.g., silver), the force spec¬ 
trum is proportional to the square of the polarization, 
i.e., F cx Re{ata*) cx jocAgP, and is absorptive when 
Re{aAg) crosses zero as for the case for the two plas¬ 
monic silver spheres. 

As shown above, the optical gradient force based on 
a collective plasmonic resonance can reach few pN. The 
force based on a single particle excitation in the form 
of a molecular resonance is weaker by three orders of 
magnitude (few fN) under similar tip-sample geometric 
conditions. In the following we address the question of 
AFM sensitivity for its possible detection and distinction 
from thermal expansion. The sensitivity of force detec¬ 
tion in AFM is limited by thermal fluctuations of the 
cantilever [3l|24]. The smallest detectable force is when 
the cantilever oscillation amplitude equals that induced 
by thermal noise given by 


AkeTkE 

ujQ 


1/2 


(3) 


with Boltzmann constant fcs, temperature T, spring con¬ 
stant A:, detection bandwidth F, and Q factor of the can¬ 
tilever. The minimal detectable force gradient due to a 
change in sample-tip interaction operating in non-contact 
mode can also be derived as [51 [53] 


F' ■ 

min 


1 hlkeTkE 
Aq y ojQ 


(4) 


with cantilever oscillation amplitude Aq. Using typical 
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FIG. 4. Simulated optical gradient force (a) and absorption (b) spectra for two silver spheres of r = 10 nm at different 
separations. The peak magnitude of optical force on resonance increases by a factor ~ 10^ from z = 20 nm to 2 = 1 nm. 
Meanwhile the resonance of the force spectrum red shifts due to plasmonic coupling. In comparison, the magnitude of the 
corresponding absorption spectra has no significant distance dependence, (c) The optical gradient force, both on-resonance 
(red) and off-resonance (blue), follows an power law distance dependence on z. 


values from force detection experiments [n na [m [THi 
[25] . we get F„iin = 0.3 pN and = 0.02 pN/nm. 

This makes collective particle excitation induced force 
detectable under typical laser intensities below the sam¬ 
ple damage threshold. In contrast force detection due 
to molecular electronic and vibrational resonance is two 
orders of magnitude below the detection limit of conven¬ 
tional room temperature AFM. 

The effect due to thermal expansion Ah can be esti¬ 
mated based on Aft, = arAT m- Assuming a sample 
thickness of 10 nm, typical thermal expansion coefhcients 
ar = 10“^ — 10“®/K, and a sample temperature increase 
of a few K after laser illumination nil US], we get a ther¬ 
mal expansion value of Aft = 0.1 - 100 pm. When this 
expansion of the sample is modulated at a resonance fre¬ 
quency of the cantilever, the corresponding force becomes 
^thermal = kAh X Q, where k is the cantilever spring con¬ 
stant of typically 3 N/m, and Q =100 - 200. The re¬ 
sulting force is in the range of ^thermal = 10 pN - 1 nN. 
Thus the thermal expansion effect can readily be mea¬ 
sured with AFM as shown experimentally even down to 
monolayer sensitivity with a single chemical bond expan¬ 
sion of few pm la¬ 
in previous attempts to detect the optical gradient 
force, a large force magnitude of F" ^ 2 pN and an ex¬ 
cellent agreement of the force spectrum with the far-held 
molecular absorption spectrum [Hill] strongly suggest 
that the observed experimental results were in fact due 
to thermal expansion and not due to the optical gradi¬ 
ent force as claimed. Optical gradient force detection has 
also been assigned to the force contrast on bowtie anten¬ 
nas reaching few pN on resonance [18) . This assignment 
is feasible in principle according to our predictions. How¬ 
ever, a force due to thermal expansion of similar magni¬ 
tude is expected, leaving the underlying mechanism of 
the imaging contrast unclear. 

In practice, the differences both in distance dependence 
and spectral frequency dependence of the optical gradi¬ 
ent force and thermal expansion can be used to differen¬ 
tiate the two mechanisms. Notably, the optical gradient 


force is a longer range effect, determined by the spatial 
extent of the near-held of the tip apex, and thus follows 
a continuous change with distance as determined by the 
tip radius (Fig. and Fig. |^). In contrast, the force 
exerted on the tip through sample expansion requires a 
direct physical tip-sample contact, and should decay on 
even sub-nm distances above the sample, independent of 
tip radius. 

A way to differentiate the two mechanisms in imag¬ 
ing plasmonic optical antennas is to map the spatial 
force distribution. The optical force is proportional to 
the local optical electric held, while the thermal ex¬ 
pansion/absorption is due to resistive heating associated 
with electric currents. Thus distinct spatial maps result 
for the two different mechanisms since the current distri¬ 
bution peaks at the positions of minima in the electric 
held in optical antennas m- 

Fundamentally, the optical gradient force due to op¬ 
tical polarization and thermal expansion due to energy 
dissipation are two dynamic processes occurring on differ¬ 
ent time scales. This provides an additional opportunity 
for their distinction. Optical polarization with coherent 
excitation is induced nearly instantaneously in the fs to 
ps range, as determined by the spectral line width. In 
contrast, thermalization of an optical excitation under¬ 
lying optical absorpton leads to thermal expansion on ns 
time scales. 

In summary, the dispersive line shape in probing elec¬ 
tronic and vibrational materials resonances as related 
to the real part of the dielectric function of the sample 
can serve as a distinguishing attribute in optical gradient 
force microscopy. Albeit weak with forces in the fN to pN 
range, their detection with advanced atomic force tech¬ 
niques can provide for a novel form of optical scanning 
probe nano-imaging and -spectroscopy. Compatible with 
a wide range of spectroscopies, including coherent and ul¬ 
trafast techniques, it can complement related all-optical 
scattering scanning near-held microscopy (s-SNOM). Our 
theory provides for the general description of the key pa¬ 
rameters of strength, distance dependence, and spectral 
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behavior for a simplified tip geometry, yet can readily 
be extended to specific tip geometries and other opti¬ 
cal processes, including inelastic, e.g., Raman, as well as 
nonlinear excitations. 


ACKNOWLEDGEMENT 

Funding from the National Science Foundation (NSF 
Grant CHE 1306398) is gratefully acknowledged. We 
thank Prof. M. A. Belkin and H. K. Wickramasinghe 
for valuable discussions. 


* markus.raschke@colorado.edu 

[1] G. Anetsberger, O. Arcizet, E. Gavartin, Q. Unterreith- 
meier, E. Weig, J. Kotthaus, and T. Kippenberg, Nature 
Physics 5, 909 (2009) 

[2] J. Berthelot, S. S. Acimovic, M. L. Juan, M. P. Kreuzer, 
J. Renger, and R. Quidant, Nature nanotechnology 9, 
295 (2014) 

[3] F. J. Giessibl, Reviews of Modern Physics 75, 949 (2003) 

[4] M. Emmrich, F. Huber, F. Pielmeier, J. Welker, 
T. Hofmann, M. Schneiderbauer, D. Meuer, S. Polesya, 
S. Mankovsky, D. Kodderitzsch, H. Ebert, and F. J. 
Giessibl, Science 348, 308 (2015) 

[5] E. Betzig, J. K. Trautman, T. D. Harris, J. S. Weiner, 
and R. L. Kostelak, Science 251, 1468 (1991) 

[6] F. Zenhausern, Y. Martin, and H. K. Wickramasinghe, 



[9] D. C. Kohlgraf-Owens, L. Greusard, S. Sukhov, Y. D. 
Wilde, and A. Dogariu, Nanotechnology 25, 035203 
(2014) 


[10] N. Satoh, T. Fukuma, K. Kobayashi, S. Watanabe, T. Fu¬ 
ji!, K. Matsushige, and H. Yamada, Applied Physics Let¬ 
ters 96, 3 (2010) 

[11] 1. Rajapaksa, K. Uenal, and H. K. Wickramasinghe, 
Applied Physics Letters 97, 2010 (2010) 

[12] 1. Rajapaksa and H. Kumar Wickramasinghe, Applied 
Physics Letters 99, 161103 (2011) 

]13] F. De Angelis, R. P. Zaccaria, and E. Di Fabrizio, Optics 
Express 20, 29626 (2012) 

]14] J. Jahng, J. Brocious, D. A. Fishman, F. Huang, X. Li, 
V. A. Tamma, H. K. Wickramasinghe, and E. O. Potma, 
Physical Review B 90, 155417 (2014) 

[15] A. Dazzi, R. Prazeres, F. Glotin, and J. M. Ortega, 
Optics Letters 30, 2388 (2005) 

[16] A. Dazzi, C. B. Prater, Q. Hu, D. B. Chase, J. F. Rabolt, 
and C. Marcott, Applied Spectroscopy 66 N62, 1365 
( 2012 ) 

[17] F. Lu, M. Jin, and M. A. Belkin, Nature Photonics 8, 
307 (2014) 

]18] J. Kohoutek, D. Dey, A. Bonakdar, R. Gelfand, A. Sklar, 
O. G. Memis, and H. Mohseni, Nano Letters 11, 3378 
( 2011 ) 

[19] A. M. Katzenmeyer, J. Chae, R. Kasica, G. Holland, 
B. Lahiri, and A. Centrone, Advanced Optical Materials 
2, 718 (2014) 

[20] J. Jahng, J. Brocious, D. A. Fishman, S. Yampolsky, 
D. Nowak, F. Huang, V. A. Apkarian, H. K. Wickramas¬ 
inghe, and E. O. Potma, Applied Physics Letters 106, 
083113 (2015) 

]21] F. Huang, V. Ananth Tamma, Z. Mardy, J. Burdett, 
and H. K. Wickramasinghe, Scientific Reports 5, 10610 
(2015) 

[22] S. Berweger, D. M. Nguyen, E. A. Muller, H. A. Bech¬ 
tel, T. T. Perkins, and M. B. Raschke, Journal of the 
American Chemical Society 135, 18292 (2013) 

[23] I. Romero, J. Aizpurua, G. W. Bryant, and F. J. Garcia 
De Abajo, Optics Express 14, 9988 (2006) 

[24] Y. Martin, C. C. Williams, and H. K. Wickramasinghe, 
Journal of Applied Physics 61, 4723 (1987) 

[25] F. Lu and M. A. Belkin, Optics Express 19, 19942 (2011) 

[26] A. Downes, D. Salter, and A. Elfick, Optics express 14, 
5216 (2006) 

[27] R. L. Olmon, M. Rang, P. M. Krenz, B. A. Lail, L. V. 
Saraf, G. D. Boreman, and M. B. Raschke, Physical 
Review Letters 105, 167403 (2010) 









